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A  Functional  for  Finite  Element 
Analysis  of  Low  Intensity 
Magnetic  Fields  Including 
Permanent  Magnetization  Effects 


Introduction 

The  analysis  of  static  magnetic  fields  using  scalar  potential  functions  in 
association  with  the  Biot-Savart  law  has  been  reported  in  the  open  literature  by 
Zienkiewicz  et  al.1  and  Armstrong  et  al.2  Wikswo3  has  indicated  how  widely 
distributed  current  distributions  can  be  economically  treated  using  this  approach. 
Although  the  work  of  Zienkiewicz  et  al.1  uses  a  general  variational,  there  is  no 
discussion  on  how  a  problem  with  permanent  magnetization  effects  should  be 
approached. 


Solution  Approach 

If  the  magnetic  fields  are  of  low  intensity  or  do  not  change  rapidly  over  the 
magnetized  body,  the  change  in  magnetic  flux  as  a  function  of  magnetic  field  in¬ 
tensity  will  be  linear;  i.e.. 


B  =  Br  +  maH,  (1) 

where  nA  is  the  incremental  permeability  and  BR  is  the  remanent  magnetization 
measured  by  drawing  a  tangent  from  the  operating  point  on  the  B-H  curve  to  the  H 
=  Oaxis. 

Following  the  usual  scalar  potential  solution  to  magnetic  Held  problems,  we 
have 


H  =  Hc  +  H,  +  f* ,  (2) 

where  Ffc  is  the  Biot-Savart  law  magnetic  field  strength  due  to  currents,  If,  is  the 
applied  farfield  magnetic  field  strength,  and  +  is  an  unknown  continuous  scalar 
potential  function.  For  the  purposes  of  this  discussion,  Ifc  and  if,  are  assumed  to 
have  either  low  amplitudes  or  to  be  changing  slowly,  so  that  can  be  taken  as  a 
constant  for  the  problem. 

Using  the  assumed  solution  of  equation  (2),  we  immediately  satisfy  one  of 
Maxwell’s  equations  involving  the  curl  of  the  magnetic  field  strength  H  and  the 
current  density T. 

?xIT-7.  (3) 
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Furthermore,  the  continuity  requirements  on  the  scalar  potential  function 
provide  for  the  satisfaction  of  the  boundary  condition  between  dissimilar  materials 
involving  the  magnetic  field  strength  components  tangent  to  the  boundary: 

ITxHI,  ="nxH|2,  (4) 

wherein  is  a  vector  normal  to  the  boundary  and  the  symbols  |,  and  |2  indicate  that 
the  operation  is  to  be  carried  out  in  materials  1  and  2,  respectively. 

The  relation  requiring  that  the  divergence  of  the  flux  density  must  vanish  fur¬ 
nished  the  field  relation  for  the  determination  of  the  scalar  potential  function: 

V*"B  =  0.  (5) 


Substitution  of  equations  (1)  and  (2)  into  equation  (5)  gives 

v- pjv*  =  -? •  M4HC  -  V  •  |i4H,  -1  •  Br  .  (6) 

The  boundary  condition  requiring  that  the  normal  component  of  the  flux  density 
be  constant  across  the  boundary  must  also  be  satisfied: 


n*  B|,  =  n«  B|2.  (7) 

A  further  boundary  condition  on  4  is  that  it  decay  to  zero  as  the  distance  from  the 
permanently  magnetized  body  becomes  large: 

4  =  0as|f-FB|-*oo,  (8) 

where  |F-7B|  is  the  distance  between  the  point  r  at  which  4  is  evaluated  and  the 
nearest  point  of  the  finitely  sized  permanently  magnetized  region. 


Functional  and  Finite  Element  Equation 

A  functional  that,  when  minimized,  leads  to  the  Add  and  boundary  conditions  of 
equations  (6)  and  (7)  is  given  below. 

X  -  f  1/2 (74  +  H,  +  Hc) •  (^4  +  mpH,  +  m4Hc  +  2B,i)dV .  (9) 

This  functional  may  be  used  as  the  basis  of  a  finite  element  formulation.  In  such  a 
formulation  the  scalar  potential  function  is  set  equal  to  the  sum  of  the  products  of 
interpolation  functions  with  discrete  unknowns: 

♦  -i"i4„  00) 
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where  N(  =  Ns  (x,  y,  z)  are  the  interpolation  functions,  are  discrete  unknown 
values  of  potential  associated  with  various  discrete  points  in  the  problem  domain, 
and  n  is  the  number  of  discrete  unknowns. 

A  discussion  of  the  details  of  the  finite  element  method  is  unnecessary  here  since 
such  details  are  covered  in  any  standard  text  on  the  subject.4 

When  equation  (10)  is  substituted  into  the  variational  principle  (9)  and  the 
functional  is  minimized  with  respect  to  the  4j  unknowns,  a  set  of  linear  algebraic 
equations  results  that  can  be  conveniently  represented  in  matrix  form: 


CK]  {4}  +  (fc>  +  {f,}  +  {fR}  =  0,  (11) 

where 

kjj  =  f"VNt*  M^NjdV  (12) 

fCi=  pN-M^dV  (13) 

1 

f,.=  *  MaHj  dV  (14) 

I  “'y 

fR  =  J^N.*BRdV.  (15) 


Equations  (12)  and  (13)  are  identical  to  those  given  by  Zienkiewicz  et  al.' 
Equation  (14)  shows  the  “loading”  term  relation  that  accounts  for  the  incident 
magnetic  field  on  the  body.  If  we  use  this  solution  approach  for  a  problem  with  an 
incident  magnetic  field,  the  finite  element  solution  gives  us  the  magnetic  field 
disturbance  caused  by  the  body  directly.  Often  the  magnetic  field  disturbance  is  the 
topic  of  interest  in  such  problems.  Equation  (15)  shows  the  term  accounting  for 
permanent  magnetization  effects;  the  distribution  of  remanent  flux  density  BR  must 
be  known  for  this  formulation. 


Conclusions 

Equations  (11)  through  (15)  allow  determination  of  the  magnetic  field  around  a 
permeable  body.  The  equations  account  for  the  effects  of  currents,  incident  fields, 
and  permanent  magnetization.  Since  the  heat  transfer  elements  in  standard  finite 
element  programs  have  the  same  formulation  as  that  given  by  equation  (12),  such 
programs  may  be  used  to  solve  a  variety  of  magnetic  field  problems.  A  preprocessor 
must  be  developed,  however,  to  calculate  the  loadings  defined  by  equations  (13) 
through  (15),  and  a  postprocessor  must  be  used  to  perform  the  additional  operation 
given  by  equation  (2)  in  order  to  obtain  the  total  field  solution. 
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